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TITLE OF THE INVENTION 

METHOD FOR PREDICTING A 
GOLFER'S BALL STRIKING PERFORMANCE 
(Corporate Docket Number PU2168) 

CROSS REFERENCES TO RELATED APPLICATIONS 
This application is a continuation application of co-pending U.S. Patent Application 
Number 10/248,332, filed on January 9, 2003, which is a continuation of U.S. Patent Application 
Number 09/683,396 filed on December 21, 2001, now U.S. Patent Number 6,506,124. 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

Not Applicable 

BACKGROUND OF THE INVENTION 
Field of the Invention 
The present invention relates to a method for predicting a golfer's ball striking 
performance for a multitude of golf clubs and golf balls. More specifically, the present invention 
relates to a method for predicting a golfer's ball striking performance for a multitude of golf 
clubs and golf balls without the golfer actually using the multitude of golf clubs and golf balls. 

Description of the Related Art 
For over twenty-five years, high speed camera technology has been used for gathering 
information on a golfer's swing. The information has varied from simple club head speed to the 
spin of the golf ball after impact with a certain golf club. Over the years, this information has 
fostered numerous improvements in golf clubs and golf balls, and assisted golfers in choosing 



golf clubs and golf balls that improve their game. Additionally, systems incorporating such high 
speed camera technology have been used in teaching golfers how to improve their swing when 
using a given golf club. 

An example of such a system is U.S. Patent Number 4,063,259 to Lynch et al, for a 
Method Of Matching Golfer With Golf Ball, Golf Club, Or Style Of Play, which was filed in 
1975. Lynch discloses a system that provides golf ball launch measurements through use of a 
shuttered camera that is activated when a club head breaks a beam of light that activates the 
flashing of a light source to provide stop action of the club head and golf ball on a camera film. 
The golf ball launch measurements retrieved by the Lynch system include initial velocity, initial 
spin velocity and launch angle. 

Another example is U.S. Patent Number 4,136,387 to Sullivan, et al., for a Golf Club 
Impact And Golf Ball Launching Monitoring System, which was filed in 1977. Sullivan 
discloses a system that not only provides golf ball launch measurements, it also provides 
measurements on the golf club. 

Yet another example is a family of patent to Gobush et al, U.S. Patent Numbers 
5,471,383 filed on September 30, 1994; 5,501,463 filed on February 24, 1994; 5,575,719 filed on 
August 1, 1995; and 5,803,823 filed on November 18, 1996. This family of patents discloses a 
system that has two cameras angled toward each other, a golf ball with reflective markers, a golf 
club with reflective markers thereon and a computer. The system allows for measurement of the 
golf club or golf ball separately, based on the plotting of points. 

Yet another example is U.S. Patent Number 6,042,483 for a Method Of Measuring 
Motion Of A Golf Ball. The patent discloses a system that uses three cameras, an optical sensor 
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means, and strobes to obtain golf club and golf ball information. 

However, these disclosures fail to provide a system or method that will predict a golfer's 
performance with a specific golf club or golf ball in different atmospheric conditions, without 
having the golfer physically strike the specific golf ball with the specific golf club. More 
specifically, if a golfer wanted to know what his ball striking performance would be like when he 
hit a CALLAWAY GOLF® RULE 35® SOFTFEEL™ golf ball with a ten degrees 
CALLAWAY GOLF® BIG BERTHA® ERC® II forged titanium driver, the prior disclosures 
would require that the golfer actually strike the CALLAWAY GOLF® RULE 35® 
SOFTFEEL™ golf ball with a ten degrees CALLAWAY GOLF® BIG BERTHA® ERC® H 
forged titanium driver. Using the prior disclosures, if the golfer wanted to compare his or her 
ball striking performance for ten, twenty or thirty drivers with one specific golf ball, then the 
golfer would have use each of the drivers at least once. This information would only apply to the 
specific golf ball that was used by the golfer to test the multitude of drivers. Now if the golfer 
wanted to find the best driver and golf ball match, the prior disclosures would require using each 
driver with each golf ball. Further, if the golfer wanted the best driver/golf ball match in a 
multitude of atmospheric conditions (e.g. hot and humid, cool and dry, sunny and windy, . . .etc.) 
the prior disclosures would require that the golfer test each driver with each golf ball under each 
specific atmospheric condition. 

« 

Thus, the prior disclosures fail to disclose a system and method that allow for predicting a 
golfer's ball striking performance for a multitude of golf clubs and golf balls without the golfer 
actually using the multitude of golf clubs and golf balls. 
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BRIEF SUMMARY OF THE INVENTION 



It is thus an object of the present invention to provide a system and method that 
allow for predicting a golfer's ball striking performance for a multitude of golf clubs and golf 
balls without the golfer actually using the multitude of golf clubs and golf balls. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 
FIG. 1 is a flow chart of the general method of the present invention. 
FIG. 1 A is a flow chart illustrating the inputs for the golf club head properties. 
FIG. IB is a flow chart illustrating the inputs for the golf ball properties. 
FIG. 1C is a flow chart illustrating the inputs for the pre-impact swing properties. 
FIG. ID is a flow chart of the inputs for the ball launch parameters. 
FIG. IE is a flow chart of the outputs that are generated for the predicted performance. 
FIG. 2 is a perspective view of the monitoring system of the present invention. 
FIG. 3 is a front view of a golf club with markers for use in determining the pre-impact 
properties. 

FIG. 3 A is a graphic of global coordinates of the markers on the golf club of FIG. 3. 
FIG. 4 is an image frame of a golfer's swing composed of a multitude of pre-impact 
exposures. 

FIG. 5 illustrates an input screen. 

FIG. 6 is an illustration of markers of a golf club on a three-dimensional plot for six pre- 
impact exposures. 

4 



FIG. 6A is a three-dimensional plot of the extrapolated head position and orientation. 

FIG. 6B is a graphic of global coordinates of the markers of FIG. 6. 

FIG. 7 is a graphic of an input menu for impact locations. 

FIG. 8 is a flow chart of the components of the pre-swing properties of FIG. 1. 

FIG. 9 is a table of the image times (in microseconds) of FIG. 8 for Golfer A and Golfer 

B. 

FIG. 10 is a table of the measured points (in millimeters) of FIG. 8 for Golfer A and 
Golfer B. 

FIG. 1 1 is a table of the static image points (in millimeters) of FIG. 8 for Golfer A and 
Golfer B. 

FIG. 12 is a table of the golf club head properties of FIGS. 1 and 1A for Golfer A and 
Golfer B. 

FIG. 13 is a table of the pre-impact swing properties of FIGS. 1 and 1C for Golfer A and 
Golfer B. 

FIG. 14 is a table of the golf ball properties of FIGS. 1 and IB for Golfer A and Golfer B. 
FIG. 15 is a table of the ball launch parameters of FIGS. 1 and ID for Golfer A and 
Golfer B. 

FIG. 16 is a table of the atmospheric conditions of FIG. 1 for a warm day and a cold day. 
FIG. 17 is a table of the predicted performance of FIGS. 1 and IE for Golfer A and Golfer 

B. 
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DETAILED DESCRIPTION OF THE INVENTION 



As shown in FIG. 1, a method for predicting a golfer's ball striking performance is 
generally designated 200. The method 200 commences with inputting information on a specific 
golf club, specific golf ball, and the swing characteristics of a golfer. At block 202, the club head 
properties of the specific golf club are selected from a database of stored and previously collected 
club head information. The specific information for the club head properties is set forth in 
greater detail below. At block 204, the pre-impact swing properties of the golfer are collected 
and stored in a database. The specific information for the golfer's pre-impact swing properties is 
set forth in greater detail below. At block 206, the golf ball properties of the specific golf ball are 
selected from a database of stored and previously collected golf ball information. The specific 
information for the golf ball properties is set forth in greater detail below. 

At block 208, the information from blocks 202, 204 and 206 are inputted into a rigid body 
code. The rigid body code is explained in greater detail below. At block 210, the rigid body 
code is used to generate a plurality of ball launch parameters. At block 212, information 
concerning the atmospheric conditions is selected from a database of stored atmospheric 
conditions. At block 214, information concerning the lift and drag properties of the golf ball are 
collected and stored. The lift and drag properties of golf balls are measured using conventional 
methods such as disclosed in U.S. Patent Number 6,186,002, entitled Method For Determining 
Coefficients Of Lift And Drag Of A Golf Ball, which is hereby incorporated by reference in its 
entirety. The lift and drag coefficients of a number of golf balls at specific Reynolds numbers are 
disclosed in U.S. Patent Number 6,224,499, entitled A Golf Ball With Multiple Sets Of Dimples, 
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which pertinent parts are hereby incorporated by reference. 

At block 216, the ball launch parameters, the atmospheric conditions and the lift and drag 
properties are inputted into a trajectory code. At block 218, the trajectory code is utilized to 
predict the performance of the golfer when swinging the specific golf club, with the specific golf 
ball under the specific atmospheric conditions. Trajectory codes are known in the industry, and 
one such code is disclosed in the afore-mentioned U.S. Patent Number 6,186,002. The USGA 
has such a trajectory code available for purchase. 

FIG. 1 A is a flow chart illustrating the inputs for the golf club head properties of block 
202. The measurements for the face properties are collected at block 401 . The face properties 
include the face geometry, the face center, the bulge radius and the roll radius. The 
measurements for the mass properties of the golf club head are collected or recalled from a 
database at block 402. The mass properties include the inertia tensor, the mass of the club head, 
and the center of gravity location. The measurement for the coefficient of restitution of the golf 
club head using a specific golf ball is collected at block 403. The measurements for the loft and 
lie angles of the golf club head are collected at block 404. The data collected at blocks 401-404 
is inputted to create the golf club head properties at block 202 of FIG. 1. 

FIG. IB is a flow chart illustrating the inputs for the golf ball properties of block 206. 
The measurement of the mass of the golf ball is collected at block 405. The measurement of the 
radius of the golf ball is collected at block 406. The measurement of the moment of inertia of the 
golf ball is collected at block 407. The measurement of the coefficient of restitution of the golf 
ball is collected at block 408. The data collected at blocks 405-408 is inputted to create the golf 
ball properties at block 206 of FIG. 1. 
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FIG. 1C is a flow chart illustrating the inputs for the pre-impact swing properties of 
block 204. The measurement of the linear velocity of the golf club being swung by the golfer is 
collected at block 409. The measurement of the angular velocity of the golf club being swung by 
the golfer is collected at block 410. The measurement of the golf club head orientation is 
collected at block 41 1. The information of the club head impact location with the golf ball is 
determined at block 412. The data collected at blocks 409-412 is inputted to create the pre- 
impact swing properties at block 204 of FIG. 1. 

FIG. ID is a flow chart of the inputs for the ball launch parameters at block 214 of FIG. 1 . 
The post impact linear velocity of the golf ball is calculated at block 416. The post impact 
angular velocity of the golf ball is calculated at block 417. The launch angle of the golf ball is 
calculated at block 418. The side angle of the golf ball is calculated at block 419. The speed of 
the golf ball is calculated at block 420. The spin of the golf ball is calculated at block 421. The 
spin axis of the golf ball is calculated at block 421. The information from blocks 416-421 is 
inputted to the ball launch parameters at block 214 of FIG. 1. 

FIG. IE is a flow chart of the outputs from the trajectory code that are generated for the 
predicted performance of block 218 of FIG. 1. Block 422 is the predicted total distance of the 
golf ball if struck with a specific golf club by a golfer. Block 423 is the predicted total dispersion 
of the golf ball if struck with a specific golf club by a golfer. Block 424 is the predicted 
trajectory shape (available in 3D or 2D) of the golf ball if struck with a specific golf club by a 
golfer. Block 425 is the predicted trajectory apex of the golf ball if struck with a specific golf 
club by a golfer. 

The golf club head properties of block 202 that are collected and stored in the system 
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include the mass of the golf club head, the face geometry, the face center location, the bulge 
radius of the face, the roll radius of the face, the loft angle of the golf club head, the lie angle of 
the golf club head, the coefficient of restitution ("COR") of the golf club head, the location of the 
center of gravity, CG, of the golf club head relative to the impact location of the face, and the 
inertia tensor of the golf club head about the CG. 

The mass, bulge and roll radii, loft and lie angles, face geometry and face center are 
determined using conventional methods well known in the golf industry. The inertia tensor is 
calculated using: the moment of inertia about the x-axis, Ixx; the moment of inertia about the y- 
axis, Iyy; the moment of inertia about the z-axis, Izz; the product of inertia Ixy; the product of 
inertia Izy; and the product of inertia Izx. The CG and the MOI of the club head are determined 
according to the teachings of co-pending U.S. Patent Application Number 09/916,374, entitled 
High Moment of Inertia Composite Golf Club, filed February 27, 2001, assigned to Callaway 
Golf Company, the assignee of the present application, and hereby incorporated by reference in 
its entirety. The products of inertia Ixy, Ixz and Izy are determined according to the teachings of 
co-pending U.S. Patent Application Number 09/916,374, entitled Large Volume Driver Head 
with High Moments of Inertia, filed July 26, 2001, assigned to Callaway Golf Company, the 
assignee of the present application, and hereby incorporated by reference in its entirety. 

The COR of the golf club head is determined using a method used by the United States 
Golf Association ("USGA") and disclosed at www.usga.org . or using the method and system 
disclosed in U.S. Patent Number 6,585,605, entitled Measurement Of The Coefficient Of 
Restitution Of A Golf Club, assigned to Callaway Golf Company, the assignee of the present 
application, and hereby incorporated by reference in its entirety. However, the COR of the golf 
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club head is predicated on the golf ball, and will vary for different types of golf balls. 

The golf ball properties of block 206 that are stored and collected include the mass of the 
golf ball (the Rules of Golf, as set forth by the USGA and the R&A, limit the mass to 45 grams 
or less), the radius of the golf ball (the Rules of Golf require a diameter of at least 1.68 inches), 
the COR of the golf ball and the MOI of the golf ball. The MOI of the golf ball may be 
determined using method well known in the industry. One such method is disclosed in U.S. 
Patent Number 5,899,822, which pertinent parts are hereby incorporated by reference. The COR 
is determined using a method such as disclosed in U.S. Patent Number 6,443,858, entitled Golf 
Ball With A High Coefficient Of Restitution, assigned to Callaway Golf Company, the assignee 
of the present application, and which pertinent parts are hereby incorporated by reference. 

The pre-impact swing properties are preferably determined using an acquisition system 
such as disclosed in U.S. Patent Number 6,431,990, entitled System And Method For Measuring 
A Golfer's Ball Striking Parameters, assigned to Callaway Golf Company, the assignee of the 
present application, and hereby incorporated by reference in its entirety. However, those skilled 
in the pertinent art will recognize that other acquisition systems may be used to determine the 
pre-impact swing properties. 

The pre-impact swing properties include golf club head orientation, golf club head 
velocity, and golf club spin. The golf club head orientation includes dynamic lie, loft and face 
angle of the golf club head. The golf club head velocity includes path of the golf club head and 
attack of the golf club head. 

The acquisition system 20 generally includes a computer 22, a camera structure 24 with a 
first camera unit 26, a second camera unit 28 and a trigger device 30, a teed golf ball 32 and a 
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golf club 33. The acquisition system 20 is designed to operate on-course, at a driving range, 
inside a retail store/showroom, or at similar facilities. 

The first camera unit 26 includes a first camera 40 and flash units 42a and 42b. The 
second camera unit 28 includes a second camera 44 and flash units 46a and 46b. A preferred 
camera is a charged coupled device ("CCD") camera available from Wintriss Engineering of 
California under the product name OPSIS1300 camera. 

The trigger device 30 includes a receiver 48 and a transmitter 60. The transmitter 60 is 
preferably mounted on the frame 34 a predetermined distance from the camera units 26 and 28. 
A preferred trigger device is a laser device that transmits a laser beam from the transmitter 60 to 
the receiver 48 and is triggered when broken by a club swung toward the teed golf ball 32. The 
teed golf ball 32 includes a golf ball 66 and a tee 68. Other trigger devices such as optical 
detectors and audible detectors may be used with the present invention. The teed golf ball 32 is a 
predetermined length from the frame 34, L u and this length is preferably 38.5 inches. However, 
those skilled in the pertinent art will recognize that the length may vary depending on the 
location and the placement of the first and second camera units 26 and 28. The transmitter 50 is 
preferably disposed from 10 inches to 14 inches from the cameras 40 and 44. The receiver 48 
and transmitter 60, and hence the laser beam, are positioned in front of the teed ball 32 such that 
a club swing will break the beam, and hence trigger the trigger device 30 prior to impact with the 
teed ball 32. As explained in greater detail below, the triggering of the trigger device 30 will 
generate a command to the first and second camera units 26 and 28 to begin taking exposures of 
the golf club 33 prior to impact with the teed golf ball 32. The data collected is sent to the 
computer 22 via a cable 62, which is connected to the receiver 48 and the first and second camera 
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units 26 and 28. The computer 22 has a monitor 64 for displaying an image frame generated by 
the exposures taken by the first and second camera units 26 and 28. The image frame is the field 
of view of the cameras 40 and 44. 

A first golf club 33 is preferably prepared for use with the system 20 to determine the pre- 
impact properties. Typically, the acquisition system 20 will take the average of ten swings from 
a single golfer to determine the pre-impact properties. These pre-impact swing properties will 
then be used to predict that particular golfer's performance with other golf clubs and golf balls 
under various atmospheric conditions without the golfer having to actually strike different golf 
balls with different golf clubs under various conditions. 

As shown in FIG. 3, the golf club 33 has a club head 50, a shaft 52, a face 54, scorelines 
56, a toe end 58 and a heel end 59. A plurality of markers are preferably placed on the golf club 
33 to highlight specific locations of the golf club 33. Only three marks are needed on the golf 
club to determine the pre-impact swing properties. A preferred embodiment is shown in FIG. 3. 
However, the acquisition system 20 is capable of using the basic features of the golf club 33 such 
as the scorelines, without the need for markers. A first marker 301 is placed on a tip end of the 
shaft 52. A second marker 302 is placed lower on the tip end of the shaft 52 than the first marker 
301 . A third marker 303 is placed on the high toe end 58 of the club head 50. A fourth marker 
304 is placed on a low toe end of the face 54. A fifth marker 305 is placed on a high toe end of 
the face 54. A sixth marker 306 is placed on a high heel end of the face 54. A seventh marker 
307 is placed on a low heel end of the face 54. An eighth marker 308 is placed in the center of 
the face 54. 

An image frame of the golf club 33 of FIG. 3 is created by the acquisition system 20 to 
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determine the location of the markers 304-308 or the scorelines relative to the markers 301-303. 
The loft, lie and face angle of the golf club are determined relative to the markers 301-303. This 
allows for the true golf club head 50 orientation to be measured from the markers 301-303. It is 
preferred that the markers 301-308 are highly reflective adhesive labels or be inherent with the 
golf club design. The markers 301-308 are preferred to be highly reflective since the cameras 40 
and 44 are programmed to search for two or three points that have a certain brightness such as 
200 out of a grey scale of 0-255. Two or more pre-impact exposures of the golf club 33 being 
swung by the golfer are acquired by the system 20. A preferred range of pre-impact exposures is 
three to nine, with six pre-impact exposures being the most preferred number. FIG. 5 illustrates 
an input screen to input the number and spacing of the exposures, the threshold level, the size of 
the points and the rigid relationship from the initial orientation screen. 

FIG. 4 is an image frame of four pre-impact exposures for a golfer swinging a golf club 
33. A first exposure 102a, a second exposure 102b, a third exposure 102c and a fourth exposure 
102d illustrate the golf club 33 prior to impact with the golf ball 66. The markers 301-303 are 
located in two dimensions, and then correlated in three dimensions. The marker 303 is correlated 
to the markers 301 and 302 on the shaft 52. The position of the face 54 and the tee ball 32 prior 
to impact our reconstructed and inputted to determine the pre-impact properties. 

FIG. 6 is an illustration of the markers 301, 302, 303 and 308 of a golf club 33 on a three- 
dimensional plot for six pre-impact exposures 102a-102f. The markers 301, 302, 303 and 308 
for each exposure 102a-102f are designated 301a, 301b, 301c, . . .etc. The global coordinates of 
the markers of FIG. 6 are illustrated in FIG. 6B. 

In the example of FIG. 6, the first exposure 102a is taken at 100 microseconds after the 
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trigger. The second exposure 102b is taken at 474.6 microseconds after the trigger. The third 
exposure 102c is taken at 849.3 microseconds after the trigger. The fourth exposure 102d is taken 
at 1223,9 microseconds after the trigger. The fifth exposure 102e is taken at 1598.6 microseconds 
after the trigger. The sixth exposure 102f is taken at 1973.2 microseconds after the trigger. 
In addition the location of the golf ball prior to impact is found. The ball location may be found 
prior to the player starting the back swing, assumed to be the same location from a previous shot, 
or found in the image. To determine the orientation of the golf club face 54 prior to impact the 
orientation of the markers discussed previously in FIG. 3 are oriented relative to the markers in 
FIG. 6. Where Ra and Ta are the rotation and translation matrix between 301a, 302a, 303a and 
301, 302, 303 and Rb and Tb are the rotation and translation matrix between 301b, 302b, 303b 
etc. 

[Point 308a] = [Point 308] * Ra + Ta. 

[Point 308b] = [Point 308] * Rb + Tb, etc. 
Using the equation, any point previously found on the golf club face 54 can be modeled 
from the measured points. From point 308f and the tee ball location, an estimate of the 
extrapolation time to impact can be made. Then, each series of points is curve fit with a second 
order curve fit and evaluated at the extrapolated time to give points 301g, 302g, and 303g of FIG 
6A. The extrapolated position data is used to calculate a new rotation and translation matrix and 
308g is located. Any feature on the face 54 can be rotated and translated to the impact position 
using this method and a vector normal to the face 54 created and located on the center of the face 
54. The initial impact location is defined as the location from the center of the tee ball 66 along 
the direction normal to the golf club face 54 and intersecting with the club head 50. The initial 
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impact location needs to be modified to correct for the amount that the ball will deform on the 
golf club face. A simple method is to correct the vertical impact location Vertical Correction = 
12.5 / 25.4 * sin(loft - attack angle). Lateral Correction = 12.5/25.4 * sin(face angle - path 
angle). More complex methods can be used to correct for the initial impact location. The 12.5 
mm is dependent on the swing speed of the club and is based on a 100 MPH swing. The slower 
the golf club head speed, the smaller the value. 308a - 308g and the image times are curve fit 
and Vx, Vy, and Vz are resolved for Rigid Body Code. 

Based on these six exposures 102a-102f, the predicted impact is at 2962.4 microseconds 
after the trigger. Based on this information, the pre-impact swing properties are calculated for 
the golfer. 

Once the pre-impact swing properties are determined (calculated), the rigid body code is 
used to predict the ball launch parameters. The rigid body code solves the impact problem using 
conservation of linear and angular momentum, which gives the complete motion of the two rigid 
bodies. The impulses are calculated using the definition of impulse, and the equations are set 
forth below. The coordinate system used for the impulse equations is set forth below. The 
impulse-momentum method does not take in account the time history of the impact event. The 
collision is described at only the instant before contact and the instant after contact. The force 
transmitted from the club head to the ball is equal and opposite to the force transmitted from the 
ball to the club head. These forces are conveniently summed up over the period of time in which 
the two objects are in contact, and they are called the linear and angular impulses. 

The present invention assumes that both the golf ball 66 and the golf club head 50 are 
unconstrained rigid bodies, even though the golf club head 50 is obviously connected to the shaft 
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52, and the ball 66 is not floating in air upon impact with the golf club head 50. For the golf club 
head 50, the assumption of an unconstrained rigid body is that the impact with the golf ball 66 
occurs within a very short time frame (microseconds), that only a small portion of the tip of the 
shaft 52 contributes to the impact. For the golf ball 66, the impulse due to friction between itself 
and the surface it is placed upon (e.g. tee, mat or ground) is very small in magnitude relative to 
the impulse due to the impact with the golf club head 50, and thus this friction is ignored in the 

■ 

calculations. 

In addition to the normal coefficient of restitution, which governs the normal component 
of velocity during the impact, there are coefficients of restitution that govern the tangential 
components of velocity. The additional coefficients of restitution are determined experimentally. 

The absolute performance numbers are defined in the global coordinate system, or the 
global frame. This coordinate system has the origin at the center of the golf ball, one axis points 
toward the intended final destination of the shot, one axis points straight up into the air, and the 
third axis is normal to both of the first two axis. The global coordinate system preferably follows 
the right hand rule. 

The coordinate system used for the analysis is referred to as the impact coordinate system, 
or the impact frame. This frame is defined relative to the global frame for complete analysis of a 
golf shot. The impact frame is determined by the surface normal at the impact location on the 
golf club head 50. The positive z-direction is defined as the normal outward from the golf club 
head 50. The plane tangent to the point of impact contains both the *-axis and the >>-axis. For 
ease of calculation, the x-axis is arbitrarily chosen to be parallel to the global ground plane, and 
thus theyz-plane is normal to the ground plane. The impact frame incorporates the loft, bulge 

16 



and roll of a club head, and also includes the net result of the golf swing. Dynamic loft, open or 
close to the face, and toe down all measured for definition of the impact frame. Motion in the 
impact frame is converted to equivalent motion in the global frame since the relationship 
between the global coordinate system and the impact coordinate system is known. The post 
impact motion of the golf ball 66 is used as inputs in the Trajectory Code, and the distance and 
deviation of the shot is calculated by the present invention. 
The symbols are defined as below: 
1 - (l 0 0) , the unit vector in the x-direction. 

j = (0 1 0) , the unit vector in the ^/-direction. 

k = (0 0 l) , the unit vector in the z-direction. 

m x , the mass of the club head. 

m 2 , the mass of the golf ball. 

, the inertia tensor of the club head. 

, the inertia tensor of the golf ball. 

r x = (a x b x c x ) , the vector from point of impact to the center of gravity of the club head. 

r 2 =(a 2 b 2 c 2 ) , the vector from point of impact to the center of gravity of the golf ball. 

r 3 =-r x +r 2 +<z 2 -b x +b 2 -c x +c 2 ) = (a 3 6 3 c 3 ), the vector from center of gravity 

of club head to the center of gravity of the golf ball. 
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v r =(v ,. v . . v, . . ), the velocity of the club head before impact. 

/ = ( v , 1 / v ^,i,/ v z.i,/ )> * e velocity of the club head after impact, 
v, = (v , v . . v . . ), the velocity of the golf ball before impact. 

1 } f \ y,\,t z,i,i / 

= ( v x,2,/ V ^,2,/ v z,2,/ )> * e velocity of the golf ball after impact. 
a.rL, cd „, , 6> y , , ), the angular velocity of the club head before impact. 
G) lf =(a) lf co , r , , ), the angular velocity of the club head after impact. 
S>^ =(co rl , , t . # 7 , , ), the angular velocity of the golf ball before impact. 
fi) 2/ =(o l2/ 2 r © Ji2f/ ), the angular velocity of the golf ball after impact. 



e e e 

xx xy xz 

p & e, 

xy ^yy ^yz 

e e e 

xz yz zz 



, the coefficient of restitution matrix. 



[l] = mv , definition of linear momentum. 
[h] = [i]® , definition of angular momentum. 



Conservation of linear momentum: 



m x v XJ +m 2 v 2J = m x v u + m 2 v 2i Bl 
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Conservation of angular momentum: 



[4<y 1 , / +[4*>2./ + " J 1 



~ c i v yXf 






-C2 V y,l,f +b 2 V :Xf 






+ m 2 


c 2 v x , 2 , f - a 2 v, 2J 




. a * v yXf ■ 


-b^xf . 




_ a 2 v y 2 >f -b 2 v x 2 f 


















+ m 2 


C 2 V *,2,i- a 2 V z,2,i 




. a i V yXi ■ 
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The definition of coefficients of restitution: 



(v,. w + r • («w * (" ^2 )))- kw + T • {&u * M ))) 
- W (v, M + J ■ fe,, x (" h )))" (v,, w + 7 • (<5u x (" 3 ))) 

k,2,, + * • K, X (" ? 2 )))" ( V «.U + * • X (" F l ))). 

frw + i ■ fas * (- f 2 )))- (v,,, j + i •(<»,,/ x (" 3 ))) 

(v,, v + J-K/ *(-r 2 )))-Gw +7-K/ x(-^)| 

.(v ii2>/ + * • (<B 2>/ x (- r 2 )))- (v iX/ + £ • x (- f x ))) 



B7-B9 



The tangential impulse on the ball causes both rotation and translation: 



C 2 (Vy,2,f - V y<2J )-b 2 (v 22 f - V 2M ) 

m, - c 2 (v x>2i/ - v xaJ )+ a 2 (v z>2 / - v,, 2>( ) 
^2(^,2,/ -v,, 2 ,,)-a 2 (v,, 2>/ -v yi2 J 









=[4 













B10-B12 



Equations B1-B12 can be combined to form a system of linear equations of the form: 



MM-M 



B13 
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where [a] , and {b} are determined from the known velocities before the impact, the mass 

properties of the golf ball 66 and golf club head 50, the impact location relative to the center of 
gravity of the golf ball 66 and the golf club head 50, and the surface normal at the point of 
impact, {x} contains all the post impact velocities (linear and angular), and is solved by pre- 

multiplying {b} by the inverse of [A\ , or any other method in solving system of equations in 

linear algebra. 

When the golf ball 66 is sitting on the tee 68, it is in equilibrium. The golf ball 66 will 
not move until a force that's greater than F m , the maximum static friction force between the golf 

ball 66 and the tee 68, is applied on the golf ball 66. 

/u s is the static coefficient of friction and g is gravity. 
For a golf ball 66 with 45 grams of mass, and a /i s of 0.3, 

F m =M s mg = (0.3X0.045X9.81) = 0.1 32 AT 

Assume this force is applied on the golf ball 66 for the duration of an impact of 0.0005 sec 
(which is an overestimation of the actual impulse), then the impulse, L, on the golf ball 66 is: 

L = (0.132X0.0005) = 0.0000662AT • s 
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This impulse, L, would cause the golf ball 66 to move at 0.00147 m/s (or 0.00483 ft/sec), and 
rotate at 8.08 rad/sec (or 77.1 rpm). Both of these numbers are small relative to the range of 
numbers normally seen for irons and woods. If the rigid body code of the present invention were 
to be applied to putters, then it would be preferable to include the friction force between the 
green and the golf ball 66 for the analysis. 



[«]■ 



e xx e xy e xz 

e xy e yy e yz 

e n e__ 

xz yz zz 



Each of the individual terms in the above matrix, e g , where i = x, y, z, mdj = x, y, z, 

relates the velocity in the /-direction to the y-direction. Each of the diagonal terms, where i=j 9 
indicate the relationship in velocity of one of the axis, x, y, or z, before and after the impact. Let 
x, y, z be the axis defined in the impact frame. The term e a includes all the energy that is lost in 
the impact in the normal direction of impact, and e yy are account for the complicated 

interaction between the golf ball 66 and the golf club head 50 in the tangential plane by 
addressing the end result. In general, the off diagonal terms e g , where are equal to zero for 

isotropic materials. 

As shown in FIG. 7, in predicting the performance of a golf ball struck by a golfer with a 
specific golf club under predetermined atmospheric conditions, an operator has the option of 
inputting an impact of the face 54 at a certain location regardless of the true location of impact. 
This allows for prediction of the performance of the golf club 33 for toe shots, heel shots and 
center shots. The type of golf ball may be selected, the type of golf club may be selected, the 
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atmospheric conditions including wind speed, direction, relative humidity, air pressure, 
temperature and the terrain may be selected by the operator to predict a golfer's performance 
using these input parameters along with the pre-impact swing properties for the golfer. 

The method of the present invention for predicting the performance of two different 
golfers, using two different golf clubs, with two different golf balls under two different 
atmospheric conditions is illustrated in FIGS. 8-17. Golfer B has a higher swing speed than 
Golfer A. Golfers A and B swing a test club 10 times for an average of the swing of each golfer. 
The predicted performances are for a golf club head 50 composed of steel and a golf club head 
composed of titanium, a 2-piece golf ball with an ionomer blend cover and a three-piece (wound) 
golf ball with a balata cover, and atmospheric conditions of a warm day and a cold day. 

FIG. 8 is a flow chart of the components of the pre-swing properties of block 204 of FIG. 
1. The components or inputs include the image times at block 203.7, the measured points at 
block 203.8 and the static imaged points at block 203.9. FIG. 9 is a table of the image times (in 
microseconds) of block 203.7 for Golfer A and Golfer B. FIG. 10 is a table of the measured 
points (in millimeters) of block 203.8 for Golfer A and Golfer B. FIG. 1 1 is a table of the static 
image points (in millimeters) of block 203.9 for Golfer A and Golfer B. 

FIG. 12 is a table of the golf club head properties of block 202 for golf club heads 50 
composed of titanium (Ti) and steel. Blocks 401-404 of FIG. 1A are included along with 
optional hosel height and Spin COR inputs. 

FIG. 13 is a table of the pre-impact swing properties of block 204 for each of the Golfers 
A and B. The table includes information for blocks 409-412 of FIG. 1C. 

FIG. 14 is a table of the golf ball properties of block 206 with information for blocks 405- 
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408 of FIG. IB. 

FIG. 15 is a table of the ball launch parameters of block 210 generated by the rigid body 
code. The table includes information for blocks 416-422 of FIG. ID. 

FIG. 16 is a table of the atmospheric conditions of block 214. 

FIG. 17 is a table of the predicted performance of block 218 which is generated by the 
trajectory code. The table includes information for blocks 422-425 of FIG. IE. 

From the foregoing it is believed that those skilled in the pertinent art will recognize the 
meritorious advancement of this invention and will readily understand that while the present 
invention has been described in association with a preferred embodiment thereof, and other 
embodiments illustrated in the accompanying drawings, numerous changes, modifications and 
substitutions of equivalents may be made therein without departing from the spirit and scope of 
this invention which is intended to be unlimited by the foregoing except as may appear in the 
following appended claims. Therefore, the embodiments of the invention in which an exclusive 
property or privilege is claimed are defined in the following appended claims. 
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